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Abstract 

In this article we study the full one-loop SUSY contributions to the lepton flavor violating Higgs decay 
h —>■ Tjl, within the context of the supersymmetric inverse sesaw model. We assume that both the 
right-handed neutrino masses, Mr, and their supersymmetric partner masses, are not far from 

the interesting Cl(TeV) energy scale, and we work with scenarios with large neutrino Yukawa couplings 
that transmit large lepton flavor violating effects. By exploring the behavior with the most relevant 
parameters, mainly Mr, and the trilinear sneutrino coupling A,^, we will look for regions of the 

parameter space where the enhancement of BR(/i —> rfZ) is large enough to reach values at the percent 
level, which could explain the excess recently reported by CMS and ATLAS at the CERN Large Hadron 
Collider. 


1 Introduction 

The discovery in 2012 of a new scalar particle at the 
LHC [1,2], whose mass has been set to m/j = 125.09 
± 0.21 (stat.) ± 0.11 (syst.) GeV [3], lays on the 
table the challenging issue of whether it is actually 
the Higgs boson from the standard model of particle 
physics (SM) or there is new physics beyond the SM 
(BSM). 

In this article, we focus on one of the new physics 
aspects of the discovered boson—the possibility of 
lepton flavor violating Higgs decays (LFVHD). In 
fact, very recently, the first direct search of the par¬ 
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ticular decay /i —>■ /rr (from now on, we will refer to 
both h ^ p,f and h —>■ r/2 decays in this shortened 
way) has been performed by the CMS Collabora¬ 
tion [4] , and an upper limit on the branching ratio of 
BR(/i — p,t) < 1.51 X 10“^ at 95% C.L. has been set. 
Additionally, CMS has also observed a slight excess 
with a significance of 2.4 standard deviations at to/j = 
125 GeV, whose best-fit branching ratio, if inter¬ 
preted as a signal, is BR{h pr) = (S.d^g ®) x 10“^. 
The ATLAS Gollaboration has just released their re¬ 
sults for the same h ^ pr decay [5] as well, focusing 
on hadronically decaying t leptons. ATLAS has re- 
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ported an upper limit of BR(/i —>■ /rr) < 1.85 x 10“^ 
at 95% C.L. in agreement with the previous CMS re¬ 
sult. Intriguingly, a small excess appears in one of the 
signal regions considered, even though it is not statis¬ 
tically significant. One way or another, the searches 
for lepton flavor violation (LFV) in the Higgs sector 
have entered into the percent level. The statistical 
significance is not enough to reach a strong conclusion 
yet, but any evidence of LFV would unquestionably 
mean a clear BSM signal due to the huge suppression 
of LFV in the SM because of the absence of flavor¬ 
changing neutral currents. 

In particular, the investigation of LFVHD is at 
present a very active field which is being studied in 
different models. LFVHD were considered for the 
first time in the context of the SM enlarged with 
three heavy Majorana neutrinos in [6] and later in 
the context of the type I seesaw model in [7], pre¬ 
dicting tiny rates due to the strong suppression from 
the large heavy right-handed neutrino masses. By 
contrast, in the context of the inverse seesaw model 
(ISS) [8] with right-handed neutrino masses at the 
O(TeV) energy scale, much larger LFVHD rates, up 
to 10“^, can be obtained [9]. In addition, LFVHD 
have been also analyzed with special attention in the 
literature within the framework of supersymmetric 
(SUSY) models [7, 10,11], finding branching ratios 
slightly larger than in the ISS case, up to 10“^. 

Here we will study the LFVHD within the context 
of the SUSY version of the ISS, which we refer to here 
as the SUSY-ISS model. In particular, we will present 
our estimate of the contribution to the BR(h —>■ r/l) 
from all the SUSY loops containing sneutrinos and 
sleptons which are typically different in the SUSY- 
ISS with respect to other SUSY models, due to the 
important effects induced by the right-handed neutri¬ 
nos and their SUSY partners with masses at ©(TeV). 
The potential increase of the LFVHD rates due to 
some of the new SUSY loops within the SUSY-ISS 
model was first pointed out and estimated in [12]. 
Other important enhancement due to SUSY loops 
have also been found in [13] for LFV lepton decay 
rates and other observables. Some phenomenological 
implications at the LHC of SUSY-ISS scenarios with 
large LFVHD rates within the same context as this 
work have been recently studied in [14]. 

In addition to performing a complete one-loop 
computation of the SUSY loops within the SUSY- 
ISS model, one of our main goals here is to analyze 
in detail if the enhancement due to the sneutrinos 
and sleptons loops can be sufficiently large as to ex¬ 
plain the LFVHD effect seen by CMS and ATLAS. 
Indeed, we will localize in this work some regions of 
the SUSY-ISS parameter space where this is possible. 


In section 2 we describe the SUSY-ISS model and in¬ 
troduce the parametrization we use to reproduce low- 
energy neutrino data. In section 3 we present the an¬ 
alytical results of our one loop calculation while we 
discuss our numerical predictions in section 4. 


2 The SUSY-ISS Model 


In this section, we briefly summarize the most rel¬ 
evant aspects for the present computation of the 
SUSY-ISS model, which is a well-known extension of 
the MSSM that can reproduce the observed neutrino 
masses and mixing. The MSSM superfield content is 
supplemented by three pairs of gauge singlet chiral 
superfields Ni and with opposite lepton numbers 
{i = 1,2,3). The SUSY-ISS model is defined by the 
following superpotential: 


W = VFmssm + EabNY^H^L^^ + NMrX + ]^XJixX , 

( 1 ) 


with £12 = 1 and 


bFMSSM =£ab 


EYeH^l} + DYdHiQ^ + UYunlQ"^ 



( 2 ) 


The generation indices have been suppressed and 
should be understood in a tensor notation as 
NYyH 2 L°‘ = Ni{Y^)ijH 2 L^■ In particular, all chi¬ 
ral superfields are left-handed, meaning that for 
D ,U ,E ,N ,X the spin 0 and spin ^ components 
are, for example in the case of E, [(br)* , (br)^]. Hi 
and H 2 are, respectively, the down-type and up-type 
Higgs bosons, defined as 



Then the soft SUSY breaking Lagrangian is given 

by 

^ ml^ + X^m\X* (4) 
-I- - v\{A^Y^)eLh^ -|- h.c. 

+ X\BxV^x)X* + VI[BrMr)X* + h.c., 


with 

_/-MSSM 
'^soft 


'^Rmle*R + ^Rmld*R + ul ml 


+ Sab{Q°‘)'^ rn^gQ^’ + SabiL‘'ymlL^ 
A 2 {XliXbXb + M2X^X^ 


+ M 3 A“A“ + h.c.) 

■eab \{u^r{AuYu)Q‘^H^2 


( 5 ) 
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+e'^j^{AeYe)L^H^ + B^H^Hl + h.c. 


During this study we will take all soft SUSY breaking 
masses to be flavor diagonal, making sure that the 
only sources of flavor violation are in the neutrino 
Yukawa coupling Y^, the lepton number conserving 
mass term Mji and the lepton number violating mass 
term Jlx- The only exception will be which re¬ 
ceives RGE-induced corrections coming from Y^, that, 
for phenomenological purposes, are given by [15] 

(Am|)„ = log > (6) 

where we take M = 10^® GeV for the rest of this 
work. 

After electroweak symmetry breaking, the neu¬ 
trino mass matrix in the basis {(vlY ^ i'r , X)'^ is 
given by 




rriD 

0 

\ 

Miss = 

toD 

0 

Mr 



1 « 

Mr 

MV 

/ 


where we have defined uid = YJ'w 2 , with V 2 =< >, 

Mr = and in order to agree with the 

definitions used in our previous article on LFV Higgs 
decays [9]. In the limit fix ^ 'mo Mr, it is possi¬ 
ble to diagonalize by blocks this matrix [16], leading 
to the 3x3 light neutrino mass matrix 

ATiight — mu Mr fixMf^^m^ , (8) 

which, in turn, is diagonalized by the PMNS matrix 
GpMNS [17]: 

^KvlNS-^lightGpMNS = , (9) 

where = diag(TOi^j , is the diagonal ma¬ 

trix that contains the masses of the three lightest 
neutrinos. Low-energy neutrino data can be repro¬ 
duced by using the following parametrization intro¬ 
duced in [9]: 

fix = Mr TTijj^ Mr. (10) 


In particular, this parametrization allows us to 
use the neutrino Yukawa couplings as input pa¬ 
rameters. As we showed in [9], the following three Y,^ 
textures 


y(i) 


0 1 
0.9 1 
1 1 


-l\ 

1 

1 ) 


( 11 ) 


y(2) 

T/i 


0 1 

1 1 

-1 1 


y(3) 


0 

-1 

0.8 


-1 

1 

0.5 


1 \ 
-1 
-1/ 
1 \ 
1 

0.5 j 


( 12 ) 


( 13 ) 


where / is a scaling factor, can lead to large t — fi 
flavor transition rates while suppressing fi — e and 
T — e flavor transition rates. We found that in the 
nonsupersymmetric ISS model, these could lead to 
large branching ratios for LFV Higgs decays, up to 
10“®, while still agreeing with other experimental 
constraints. As mentioned in the introduction, previ¬ 
ous studies have demonstrated that supersymmetric 
contributions usually enhance the LFV rates. In par¬ 
ticular, in the present SUSY-ISS model, since we con¬ 
sider a seesaw scale Mr not far from the electroweak 
scale, this low value will enhance the flavor slepton 
mixing due to the RGF-induced radiative effects by 
the large neutrino Yukawa couplings, and this mixing 
will in turn generate via the slepton loops an enhance¬ 
ment in the LFVHD rates. On the other hand, new 
relevant couplings appear, like A,^, which for right- 
handed sneutrinos with 0(1 TeV) masses may lead 
to new loop contributions to LFVHD that could even 
dominate [12]. In light of the recent GMS and AT¬ 
LAS searches [4, 5] for h ^ fir, this calls for a new 
and complete evaluation of the SUSY contributions 
to this observable in the SUSY-ISS model. 


3 Analytical results 

In this work, we perform a full one-loop diagrammatic 
computation of all relevant supersymmetric loops 
within the SUSY-ISS model for BR(iLa; —)• ik^m), 
where here and from now on refers to the three 
neutral MSSM Higgs bosons, Hj. = {h, H, A). This is 
in contrast to the previous estimate in [12] where an 
effective Lagrangian description of the Higgs medi¬ 
ated contributions to LFV processes was used, which 
was appropriate to capture the relevant contributions 
at large tan/?, and where the mass insertion approach 
was used to incorporate easily the flavor slepton mix¬ 
ing (Am|)ij, working in the electroweak basis. How¬ 
ever, an expansion up to the first order in the mass in¬ 
sertion approximation may not be appropriate for the 
type of scenarios studied here, due to the large flavor- 
nondiagonal matrix entries considered in this work. 
On the other hand, we are interested also in small 
and moderate tan /3 values, not just in the large tan /3 
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regime, and we also wish to explore more generic soft 
masses for the SUSY particles and scan over the rele¬ 
vant neutrino/sneutrino parameters, mainly Mu, A^, 
and , not focusing only on scenarios with univer¬ 
sal or partially universal soft parameters nor fixing 
the relevant parameters to one value as in [12], Thus, 
our calculation is performed instead in the mass ba¬ 
sis for all the SUSY particles involved in the loops, 
i.e, the charged sleptons, sneutrinos, charginos, and 
neutralinos. 

Before moving to the calculation, let us introduce 
the relevant interaction terms from the Lagrangian 
for the study of the LFV Higgs decays. Following 
the notation in [7], these terms are given in the mass 
basis by 


cancellation of divergences appearing in the loop con¬ 
tributions in both form factors and Fn^x in 

eq.( 16),giving, as expected, a finite result without the 
need of the renormalization procedure. Notice that 
this cancellation is not trivial and is, therefore, a good 
test of our results of the form factors in Appendix B. 
The parametrization of the LFVHD widths in terms 
of form factors remains unchanged and is given by 

T{Hx ^ iJm) = ^ - {-Ami^mi^Re{FL,xF*Rx) 




V 



mik - 


A.J a 


-gi 

-gi 

-iH, 


A 




Laj^R 


w 


b\1Pl + BZPr 




Xj Va + h.c ., 
xlia + h.c. , 




“ gBxXi 




— ^H Y° 

2 ^^xAa 


+ w]^,Pr 
d^^Ipl + d^SIpr 


X, 


x^ 


Bh,,u = —gHxi 


q(^) p _i_ c(^) p 


(14) 


where the coupling factors have been expressed in 
terms of the SUSY-ISS model parameters and are 
collected in Appendix A. 

We take into account the full set of 1-loop SUSY 
diagrams shown in figure 1. It is interesting to no¬ 
tice that, since we work in the mass basis, the set 
of diagrams contributing to LFV Higgs decays (four 
diagrams with charginos and sneutrinos in the loops, 
and four more with neutralinos and charged sleptons) 
is the same as in the SUSY type I seesaw model which 
was considered in [7]. We keep their definition of the 
form factors 


iFx = -iguek{-P2){FL,xPL + FR^xPR)ve^{p3 ), (15) 

where F^ is the decay amplitude for Fl^ —>■ ikim with 
again FI^ = {h, H, A) and pi = ps — P 2 is the in¬ 
going Higgs boson momentum. The contributions of 
the SUSY diagrams are summed in Fr ^ and Fr ^ 
according to 

FL,x=j2PL% FR,x = j2Pnl- ( 16 ) 

1=1 i=l 


where g is the SU(2)i coupling constant, is the 
mass of the Higgs boson while and are the 
masses of final state leptons. 

4 Numerical Results 

In this section we show the numerical results of the 
LFV decay rates of the lightest neutral Higgs boson, 
BR(/i —>■ t/2), as a function of the most relevant pa¬ 
rameters of the SUSY-ISS model for the full SUSY 
contribution to LFVHD, namely, Mr, A^, and 
It should be noted that, in the absence of CP vio¬ 
lation, as in our case, BR(/i —>• Tfi)=BR{h —>• g,f) 
and, therefore, in comparing with data the two rates 
should be added. 

We have imposed various experimental con¬ 
straints, choosing as example two benchmark points 
leading to a Higgs boson mass within Icr of the cen¬ 
tral value of the latest CMS and ATLAS combination, 
and with supersymmetric spectrum allowed by AT¬ 
LAS and CMS searches. Indeed, we concentrate for 
this work on the slepton sector, since the squark sec¬ 
tor is irrelevant for the LFVHD. The squark param¬ 
eters are only relevant for the Higgs mass prediction, 
thus one can always adjust the squark masses and 
the trilinear couplings A^ and Af, in order to ensure a 
correct Higgs boson mass. We also restrict ourselves 
to the case of Mr > rrih, avoiding constraints from 
the invisible Higgs decay widths, and consider only 
real Upmns and mass matrices, making constraints 
from lepton electric dipole moments irrelevant. Fi¬ 
nally, we also take into account the LFV radiative 
decays whose current upper limits at the 90% C.L. 
are 


Their analytic expressions are taken from [7] and re¬ 
produced in Appendix B for completeness, includ¬ 
ing the proper modifications to adapt them to the 
SUSY-ISS model. We have checked analytically the 


BR{^l ex) < 5.7 x 10 "^^ 

[ 18 ], 

( 18 ) 

BR(t ^ ex) < 3.3 X 10 "® 

[ 19 ], 

( 19 ) 

BR(r —)■ fix) ^ 4.4 x 10 “® 

[ 19 ]. 

( 20 ) 
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Figure 1: One-loop supersymmetric diagrams contributing to the process —>■ tk^m- 


Points excluded by LFV radiative decays will be de¬ 
noted by a cross, while a triangle will represent the 
ones allowed. We present here the predictions of 
BR(/i — >■ tJJl) for the three neutrino Yukawa tex¬ 
tures exposed in the section 2, ensuring the practi¬ 
cally vanishing LFV in the fi — e sector, i.e., leading 
to BR(/x —>• ey) ~ 0 and BR(h —)• ep,) ~ 0. It should 
be noticed that these textures also suppress substan¬ 
tially the LFV in the r—e sector. Therefore, the most 
stringent constraint, making use of these textures, is 
that of the related LFV radiative decay r —>■ //y. 

In hgure 2, we show the behavior of BR(/i —r/2) 
as a function of Mr for the three textures presented 
in the previous section, (upper left panel), Yr'f} 

/o'! 

(upper right panel), and Y^^' (lower left panel), for 
different values of the scaling factor / = O.OI, 0.1, 
I, -s/dTr, First of all, we clearly see that, as 

expected, the larger the value of / is, the larger the 
LFV rates are. We also observe qualitatively differ¬ 
ent behaviors of the LFV rates between small (/ < 1) 
and large (/ > I) neutrino Yukawa couplings. As we 
have checked, this difference comes from the different 
behavior with Mji of the two participating types of 
loops, the ones with charged sleptons where the LFV 
is generated exclusively by the mixing (Am?)ij and 
the ones with sneutrinos where the LFV is generated 
by both (Am|)ij and {Y^)ij. In the case of small /, 
charged slepton-neutralino loops dominate and they 
only depend logarithmically on as can be seen 
from eq. (6), leading to the apparent flat behavior. 
However, we checked that this flat behavior disap¬ 
pears when both Mr and Mg (and as a consequence 
all slepton and sneutrino masses) increase simulta¬ 
neously. When the scale factor / becomes larger, 
contributions from sneutrino-chargino loops become 


sizable and even dominate at low M^. They decrease 
with Mfi, due to the increase in the singlet sneutrino 
masses, which explains the decrease in BR(/i —>■ rp) 
observed in the upper plots and on the left-hand side 
of the bottom plots for large / > I. In the latter, 
the appearance of dips due to negative interferences 
between the two types of loops marks the transition 
between the two regimes, with the main contribu¬ 
tion coming from sneutrino-chargino loops at low Mji 
and from slepton-neutralino loops at large Mr. For 
the first benchmark point, the largest BR(h —>■ rp), 
allowed by the t ^ fi'y upper limit, are obtained 
for / = i/dTr or i/bTr and M/j < 2 TeV, with a 
value of 0(10“^) for the three textures, which could 
be probed in future runs of the LHC. Up to now, 
the trilinear neutrino coupling Ai, had been set to 
zero, whilst on the lower right panel of figure 2 we 
have chosen = 2.5 TeV and show the behavior of 
BR(/i — >■ Tp) with Mfi for the three textures with a 
scaling factor / = VStt. This value of A^ leads to a 
suppression of the r —> /xy decay rates while simulta¬ 
neously enhancing BR(h — >■ rp). As a consequence, 

very large LFVHD branching ratios can be obtained 
(3) 

for YrVi with low Mr ~ I TeV, allowed by t —>• /xy, 
achieving values up to 7 x 10“^. These large rates are 
very close to the percent level and within the sensi¬ 
tivity of the present experiments. 

We next study the behavior of BR(h —>■ rp) as 
function of the SUSY mass scales in a simplified sce¬ 
nario where all the SUSY masses are equal to a com¬ 
mon parameter tosusy, namely, 

wsuSY = mi = rrie = = m^ = Mq = Mi = M2. 

( 21 ) 

Fig. 3 left shows the expected decoupling behavior 
where BR(/i —?> rp) decreases when increasing the 
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Figure 2: BR{h —>• t / 2 ) as a function of for (upper left panel), (upper right panel), and Yrf,^ 
(lower left panel), with = mg = = 1 TeV, M 2 = 750 GeV, /r = 2 TeV, = 0, tan/3 = 5 

and different values of the scaling factor / = 0.01, 0.1, 1, -s/iTr, -s/Ott. On the lower right panel, the behavior 
of BR(/i —>■ r/2) as a function of is shown for the three textures with = mg = m^^ = m^^ = 1 TeV, 
M 2 = = 500 GeV, Ai, = 2.5 TeV, tan /3 = 10 and / = On all the panels, niA = 800 GeV and Mg = 

1 TeV. We set Aq = Ae = Bx = Br = 0 and the GUT inspired relation Mi =5/3 M 2 tan^ 9w in these and 
all the figures of the paper. Crosses (triangles) represent points in the SUSY-ISS parameter space excluded 
(allowed) by the r —>• /ry upper limit, BR(r —>• /ry) < 4.4 x 10“® [19]. 


heavy sparticle masses. This plot is for the particu¬ 
lar input Yr ]}, but similar behaviors (not shown) are 
obtained for the other two studied textures and 
Yrii . In this figure we have included the full predic¬ 
tions for BR(/i —>■ tJx), as well as the separated contri¬ 
butions coming only from chargino-sneutrino loops, 
i.e, diagrams (l)-(4) in fig.l, and from neutralino- 
slepton loops, i.e., diagrams (5)-(8) in fig.l. We see 
that not only the full prediction but also the sepa¬ 
rated contributions from these two subsets decrease 
with msuSY, showing that the decoupling occurs in 
both, the charginos-sneutrinos and the neutralinos- 
sleptons sectors, as expected from the decoupling 
theorem. We also see that, in this heavy sparti- 
cles scenario, the contributions from the charginos- 
sneutrinos sector dominate by many orders of mag¬ 


nitude over the ones from the neutralinos-sleptons 
sector. In order to better understand the contri¬ 
butions from the charginos-sneutrinos sector, which 
are the ones containing the new sparticles with re¬ 
spect to the MSSM, we consider next the simple 
case of Ami,, = 0, where the contributions from 
the neutralinos-sleptons sector vanish, and only the 
contributions from charginos-sneutrinos remain. We 
show in fig. 3 right the separated contributions from 
each diagram (1), (2), (3), and (4) of fig.l, and the full 
result. The contribution from diagram (1) is clearly 
subleading by several orders of magnitude and the 
contributions from the vertex correction, diagram (2), 
and the self-energies, diagrams (3) and (4), clearly 
compete in size. We also see that their interference 
is destructive, such that the full result, that decou- 
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Figure 3: BR(/i —> Tfj.) as function of the common SUSY mass parameter uigusY defined in eq.(21) for the 
Yukawa coupling matrix Yr}} with Mji = 1 TeV, / = ttia = 800 GeV, fi = 2 TeV, tan /? = 10 and 

Ai, = 2.5 TeV. Left panel: Contributions from chargino-sneutrino loops, denoted by v-x~, neutralino-slepton 
loops, denoted by and full results for BR(/i —)■ rp,). Right panel: Individual contributions from each 
chargino-sneutrino diagram (1), (2), (3), and (4) in hg.l and full result in the case of Ami,_ = 0, where the 
neutralino-slepton contributions vanish. 


pies with mguSY, manifests that a strong cancellation 
among self-energies and vertex corrections is happen¬ 
ing, as expected. Notice also that diagrams (3) and 
(4) do not decouple individually with m-suSY, but 
they do decouple when adding all the diagrams, as 
expected. 

Regarding the relevance for the searched enhance¬ 
ment in the SUSY contributions with respect to the 
trilinear coupling A^, we have found that the LFVHD 
rates are indeed very sensitive to the particular value 
of A^. Thus, we study in hgure 4 the behavior of 
BR(/i — >■ Tp) with this parameter for the two scenar¬ 
ios considered previously, with Mr = 1 TeV and the 
scaling factor / = On both plots we confirm the 

strong dependence of the LFVHD branching ratios 
with presenting deep dips in different positions 
that depend mainly on the values of YJy, /r, tua and 
tan j3. In particular, the — vl — vr coupling and the 

— vr mixing are controlled by these parameters, 
which would lead to the appearance of dips in the 
regime where contributions from sneutrino-chargino 
loops dominate. This is the case of hgure 4 and it 
is interesting to note that, for this choice of param¬ 
eters, practically all the parameter space is excluded 
by T —)■ /xy except the points within the dips and 
surrounding them, where the LFV radiative decay 
T —^ /xy suffers also a strong reduction. An inter¬ 
esting feature we found is that the location of the 
dips in BR{h — >• rp) and BR(t — > /xy) usually do not 
coincide, therefore allowing for large LFV Higgs de¬ 
cays rates, not excluded by t —>■ /xy, above 10“^ and 
within the reach of the LHC experiments. 


Finally, the dependence of the LFVHD rates on 
the new sneutrino soft SUSY breaking scalar masses, 
and is depicted in hgure 5 where these pa¬ 
rameters are varied independently from the SUSY 
scale. As when varying Mr, increasing and 
makes the singlet sneutrinos heavier and decreases 

the size of the chargino contribution. For and 
f2') 

Yrfj, which are dominated by this contribution, the 
BR(/i —>■ Tp) exhibits a strong decrease between 200 
GeV and 14 TeV, by more than hve orders of mag- 
nitude in the case of Yr^ . For Y-Vt a dip can be 
observed, due again to a cancellation between the 
chargino and neutralino contributions, with the lat¬ 
ter dominating at large For the hrst bench¬ 

mark point, the largest h ^ rp rates allowed by the 
r —>■ /xy upper limit are obtained for Y^^' with = 
200 GeV, with a maximum value of ^ 3 x 10“^, just 
one order of magnitude below the present LHC sen¬ 
sitivity. If we move our attention to the vicinity of 
the region of low values of for the second bench¬ 
mark point, we found large LFVHD rates, as dis¬ 
played on the right panel of figure 5, with Mr = 200 
GeV and A^, = 2.5 TeV. We observe a huge increase 
in BR(/i —>■ Tp) for the three Yukawa textures 
Yt‘^\ and with maximum values of ~ 4 x 10“^, 
~ 8x 10“^, and ~ 1.5x 10“^, respectively, due mainly 
to the low values of and Mr. Unfortunately, 
all the parameter space for Yr}} and Yrfi} cases is 
excluded by the r —>■ /xy upper limit. By contrast, 
most of the points for the YrV texture are in agree¬ 
ment with this upper bound, because they are located 
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Figure 4: Dependence of BR{h —)• rp.) on Aj, for the three neutrino Yukawa couplings Yrl\ Y^'f}, and Yrf}, 
with M 2 = 750 GeV, tan j3 = 5 and fi = 2 TeV (left panel) or with tan /3 = 10 and M 2 = fi = 500 GeV 
(right panel). On both panels, ttia = 800 GeV, Mg = 1 TeV, Mr = mi = mg = = 1 TeV, and 

the scaling factor / = v^Ott. Crosses (triangles) represent points in the SUSY-ISS parameter space excluded 
(allowed) by the r —)■ /ry upper limit, BR(r —>■ /iy) < 4.4 x 10“® [19]. 



Figure 5: Dependence of BR(h — t^l) on three neutrino Yukawa couplings Yr}}, Yrp\ 

and Yrfi\ with Mji = = rug = 1 TeV, M 2 = 750 GeV, /r = 2 TeV, tan/3 = 5 and A^, = 0 (left panel) or 

with Mji = 200 GeV, mi = mg = 1 TeV, M 2 = M = 500 GeV, tan,3 = 10 and = 2.5 TeV (right panel). 
On both panels, = 800 GeV, Mg = 1 TeV, and / = -s/Ott. Crosses (triangles) represent points in the 
SUSY-ISS parameter space excluded (allowed) by the r —>• /ry upper limit, BR(r —>• ^ 7 ) < 4.4 x 10“® [19]. 


in a region where the t — > /ry rates suffer a strong 
suppression as a consequence of the value set for 
in this case, A^, = 2.5 TeV. This fact allows us to 
obtain a maximum value of BR(/i —>■ rp) ^ 1.1%, 
completely within the reach of the current LHC ex¬ 
periments and large enough to explain the CMS and 
ATLAS excesses if confirmed by other experiments 
and/or future data. 


5 Conclusions 

In this article, we have presented the results of an 
updated and full one-loop calculation of the SUSY 
contributions to lepton flavor violating Higgs decays 
in the SUSY-ISS model. We found much larger con¬ 
tributions than in the type I seesaw due to the lower 
values of Mr ~ 0(1 TeV), an increased RGE-induced 
slepton mixing, and the presence of right-handed 
sneutrinos at the TeV scale, where both sleptons and 
sneutrinos large couplings transmit sizable LFV due 
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to the large Fj^/(47r) ~ 0{1) considered here. We 
showed that the branching ratio oi h ^ rp, exhibits 
different behaviors as a function of the seesaw and 
SUSY scale if it is dominated by chargino or neu- 
tralino loops. Moreover, a nonzero trilinear coupling 
Ai, leads to increased LFVHD rates. Choosing differ¬ 
ent benchmark points, we found that BR(ft, —)■ t^) of 
the order of 10“^ can be reached while agreeing with 
the experimental limits on radiative decays, provid¬ 
ing a possible explanation of the CMS and ATLAS 
excesses. While out of the scope of this work, a com¬ 
plete study including nonsupersymmetric contribu¬ 
tions in the SUSY-ISS model and a detailed analysis 
of experimental constraints beyond radiative LFV de¬ 
cays will be presented in a future article. 
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Appendices 

A Mass matrices and couplings in the SUSY-ISS model 

We present in this appendix the mass matrices and coupling factors that are relevant to our calculation of 
the LFV Higgs decays. The sneutrino mass matrix M? is defined by 

VL 


vr 
X* 

X 

where vr and X are vectors made of weak eigenstates and defined in a similar fashion, e.g. vr = 
. The 18 x 18 sneutrino mass matrix is expressed in terms of 3 x 3 submatrices, giving 

MIl 0 0 moM^ 0 

0 {Ml^Y 0 0 m*^MR, 

0 {Ml^Y 0 MR^Ji\ {BrMY)* 

0 0 BrM*^ 

0 lixMl 2{Bxfi*xV 

0 {BrMYY" Yx^h ‘^{BxAx) 

with 

= todtoIj-I-to|-I- cos2/3, (24) 
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tan/3”*” 

(25) 

^RR = + MrM]^ + , 

(26) 

Mxx = ^R^R + 9-x9*x + ; 

(27) 


where we have used the fact that is a symmetric matrix. Then, the sneutrino mass matrix is diagonalized 
using 

U^Mjil = M?= diag(m?^ ,..., mlj , (28) 

which corresponds to 


i>R 

X* 

V ^ J 


( n, \ 


= U 


(29) 


\ flis J 


The basis in eq. (22) uses the sneutrino electroweak eigenstates and their complex conjugate states, and they 
fulfill 



(30) 


(31) 

= Ul^n, , 

(32) 


since the physical sneutrinos are real scalar fields. While both eqs. (31) and (32) are equally valid, we choose 
eq. (31). The mass matrices of the other SUSY particles, namely the charginos, neutralinos, and charged 
sleptons, are the same as in the SUSY type I seesaw studied in [7] and we will use their definitions of the 
corresponding rotation matrices, which in turns were based on the conventions of [20] for the charginos and 
neutralinos. Concretely, U and V will be the matrices that rotate the chargino states and N the one that 
rotates the neutralino states. In addition, combinations of rotation matrices for the neutralinos are defined 
as 


iV' 1 = Nai cos 9w + Na 2 sin Ow , 
iV '2 = -Nai sin 9w + Na 2 cos 9w . 

As for the charged sleptons, they are diagonalized by 


(33) 


(34) 


where I = (cl, cr, fLL,flR, tr , tr)’^ are the weak eigenstates and i = {£i ,..., t'e)^ are the mass eigenstates. 

When compared with the SUSY type I seesaw, only the coupling factors and are modified. 

In the SUSY inverse seesaw, they are defined in the mass basis with diagonal charged leptons by 


A 


(e,M,r) 

Raj 


=u. 


(l,2,3)arji 


Vn- 


r'rD(i,2,3)fe 

•\/2mw sin f3 


U, 


/c+9,aU/2 ; 


9H^ 




— W {9LL.u)ikUi^Ukl3 + {9RR,v)ikU^^gXk+9,P 
+ i9i%UUlaUk+9,P + {9^%):kU*k+9,aU^,p 
+ i9i%MU*aUk+12,p + {9‘'LlXkUt+^2.aUr,p 
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(.9 


i^) \ _ _ 

LL,v)ik — 


'^Z (x) r (x 


2 COS 9w 


Sik + 


niw sin /3 


/ (x) \ {kkl j^TTl £)') (^x) 

\9RR.iy)ik — ^ —O', 


(ff 


mw sin /3 ® ’ 

(x) N, _ ('OlD^t)ik (x) I 

LR,i/)ik —t: ^ zO-) + 


2mw sin /3 




2mvK sin /? 




^ _{x) 


( 35 ) 


which are summed over the internal indices, with i,k = 1,...,3. We reproduced below the unmodified 
coupling factors from [7] (correcting a typo in for completeness in the mass basis with diagonal 

charged leptons 


4 ( 6 ,M,t) 
^Laj 


‘Ole,n,T fj 


B 


(e,/i,T) 

Loca 


B 


(e,/i,T) 

Raa 


=V2 

=V2 


me 


,N*R 




2mwcosf3 ( 1 . 3 , 5 ) 

1 


■a ]\t'* AT ^ 

smOwNai -^ 


cos9 


w 


R 


— sm9\YNai ~ 




R 


W 1 pr- nv; 

(1.3.5)a+ 2mw cos/? (2.4.6) 


R) 

(2,4,6)ct 

m, 


w 


=4 (^Ru;,v,i + 4->u;,v,i). 


=Tr2 i-"' 

^_L 

2 cos 9w L 




0 ’ 


i2 Vjl -r (J 2 yj2^ 

{sin9wN^i - cos9wN^2)io["''’N*3 


■o'f^N* 


+ {sm9wN:, - cos 9wK^){a[^^N*, + ] 


n(“) _n(“4 

^Rab —^Lah ’ 


Q\XJ _ 

‘^L,e — 


me 


cR) _ Q 


2myy cos /? ” ^ 

x)=^ 


{x)* 


9hJJb - *5 


(x) p*{^) p(^) _i_ I p^v^? j_ p''^-' 

yLL^e^la ^1/3 -^2^ yi?L,e'^2a 

I P*('^) p(^) _1_ P*(^) pC^) _1_ P*(^) p(-^) _1_ P*(^) p(^) 

'^yLL,iR-^3a ^3/3 -^4/3 yLi?,/x'^3Q' ^^4/3 yi?L,M'^4a 4 X 3 ^ 

I p*(^) p(^) _i_ p*(-^) p(^) _i_ p*(-^) p(^) _i_ p*(^) pC-^) 

yLL.T^boc -^5^ yRR,r^6ct -^6/3 yLR,r^3a yRL,r^6a ^5/3 




o*(-^) p(^) 


(rc) 


?*(^) p(-^) 


(a:) 7p*(^) p(^) 


(x) _ 


TOZ (x) 

-cr/ 


i — sin^ 0^12 ) + 


to| 


da:) 


mvv cos/? 


(£c) _ 


mz 
cos 0^/ 


,(a:) 


sm 


9w) 


mj 




(x) 

9LR,i 


(-'"1 


= -ar>Ae- 


R) ) 
^^5 M 


mw cos /? 
me 


2mw cos .d ’ 


(a?) (x)* 

9RL,e —9LR,e > 


(36) 


with 




( . 

sma 

— cos a 
y i sin .0 j 


da:) 


( COS a ^ 

sin a 
—icos /? y 


da:) 


sin(a + /?) ^ 

— cos(q; + /?) 

0 / 




— sma 
cos a 
0 
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(37) 


^ cos a'' 


4“’ = 


Sin a 


cos /3y 


(^) _ 
ag — 


(cos q;^ 
sin a 

V 0 ) 


1^) 


sin a 



— cos a 

, for iJj; = 

ijo 

—zsin /5y 




B Form factors in the SUSY-ISS model 

We present here the form factors that correspond to the diagrams of figure 1. The original calculation in 
the SUSY type I seesaw was carried by some of the authors in the mass basis and in the Feynman-’t Hooft 
gauge [7]. The only changes required to adapt the original form factors to the SUSY-ISS model are the 
sum over sneutrinos that has to be extended to the 18 mass eigenstates and the new couplings defined in 
Appendix A. In the following formulas, summation over all indices corresponding to internal propagators is 
understood. These would be a , ,5 = 1,, 18 for the sneutrinos, i, j = 1,2 for the charginos, a , /3 = 1,..., 6 
for the charged sleptons and a,b= 1 ,..., 4 for the neutralinos. 


F, 


( 1 ) 
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F-, 


F, 


F, 


F, 


( 2 ) 

L^x 

(3) 
L,x 

(4) 
L,x 

(5) 

L,x 


and 


(Bq + + 'rnl^Ci2 + rnj^^Cn — C'12)) 


IGtt^ - 

+me^mi^ {Cii + Co) (Cn - C'12-I-Co) Kl 4 

(Cii - C 12 ) {C 12 + Co) K^f -t to--to--C'o 


W9h^ 




IGtt^ 


-mc^iCii - C12) - me^Ci2 + '^j.-Co cf 


S 


(x) 

L,t„ 


2 2 


S 


» 

L,tk 


("^4) + "^4 ("^4) + ("*4 Si {mj^) + mi^ {mj ^)) 


("^D + ("iL) + "^4 {m}J + rrii^F^l^ i^ej) 
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^Bo + mqJCo + + '^\{Cii — Ui2)^ 


x, X 
^L1 


IGtt^ L 

{Cii + Co) (Cn - C'12 + Cq) 171^0(712 

-|-TOi!j^myo (Cii - C12) -I- mi^m^o (C12 -b Co) + myom^oCo 

-rrii^ (Cii - C12) - me^Ci2 ^£2 + i-lo 

"^4 i^k ) + "^4 ^Rs (^4 ) + (”^4 i^k ) + "^4 i'^k ) 

ijne^^R {mjj + 


H ^ 

II 

'''99Hja,ip 

IGtt^ 

II 

q{x) 

^L,irn. 

^k ~ "^L 

pis) _ 

^ L,x 

q{x) 

"^L “ ^k 

where. 



Bo = 


Bo , mi ., TO^_) in j , 
Bo {mjj , m?o, TO?o) in F^^'’ 


Co,11,12 — ■* 


C, 


o,ii,i 2 (m|^,m|^^,TO?^,m 2 _,TO 2 _) in j , 


C, 


0 , 11,12 


i'^k- 

2 

, mir 
J^X 

2 

, ml_ 

Xi 



i'^k- 

1 Tflir 

' J^x 


,m|o, 

.m?o) in 

^(5) 



’""Is 


, m~ ) in 
h 

p(6) 

^L^x ■ 
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The couplings and self-energies from the neutralino contributions to the form factors were defined as 


^L1 


X, x*" 

^L3 

X, 


^L6 


^L7 


~ ^Lcca^Rab^Rocb 

_ oC^fc) oC^m)* 

^ Rota^ Lab^ Loib 

_ oC-^fc) n(^) IDi^rn)* 

— ^ Raa^ Lab^ Rotb 

_ oC-^fc) n(^) oC-^m)* 

~ ^Laa^Rab^Loib 

_ r){^) oC^m)* 

~ ^Raa^Rab-^Rab 

_ r>(^fc) 7^(^) oC-^m)* 

~ ^Laa^Lab-^Lab 

_ r>(^fc) 7^(^) dC-^jti)* 

~ ^Laa-^Lab^Rab 



_ f>{^k) 

^ Raa^ Rj3a 

'Ll 


_ rji^k) 

— ^Laa^LjBa 

'L2 

a:,X° 

_ r>C^fc) r>(^rn.)* 

Laa^ R(3a 

'L3 






. o, O^TO^O , o n n 


\ TDi^k) p(-^m)* 

d'^i?aa'^i?aa ’ 

\ D(^fc) rjC^m.)* 
J^Laa^Raa ' 


(38) 


The couplings and self-energies from the chargino contributions to the form factors, k®’ ^ ^ , and 

can be obtained from the previous expressions ^ , 6^’ ^ and by using the following replacement 
rules m^o —)■ m, mi —)■ —> A^'-^, —)• , a —>■ i, and b ^ j. 

The form factors F^^,i = 1,..., 8 can be obtained from i = 1,..., 8 through the exchange L o i? in 
all places. 
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